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Abstract 

We compute the large N effective action of the 0{N + 1)/0{N) non-hnear sigma model 
CN ■ including the effect of the pion mass to order m^/f^. This action is more complex than 

r^ , the one corresponding to the chiral limit not only because of the pion propagators but also 

f^ I because chiral symmetry produce new interactions proportional to m^//^. We renormalize 

"j" ' the action by including the appropriate counter terms and find the renormalization group 

equations for the corresponding couplings. Then we estudy the unitarity propierties of the 
i-^ I scattering amplitudes. Finally our results are applied to the particular case of the linear 

I ' sigma model and also are used to fit the pion scattering phase shifts. 
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1 Introduction 

In recent years a lot of effort has been devoted to the so called Chiral Perturbation Theory 
(xPT) |l|, 0. This formalism provides an usefuU tool for the phenomenological description 
of the low-energy hadron dynamics in terms of some parameters that can be fitted from the 
experimental data. Even it is possible that some day they will be computed from the underlying 
theory which is expected to be Quantum Cromodynamics (QCD). 

As xPT can be applied without a precise formulation of this underlying theory, it has been 
also used for the parametrization of the unknown Symmetry Breaking Sector of the Standard 
Model giving rise to a model independent description of this sector that can be quite usefuU 
for the analysis of the future Large Hadron Collider (LHC) data. 

However, xPT still has some problems that should be solved in order to be completely usefuU 
for the practical applications. Usually xPT computations are done to the one-loop level, or what 
it is the same, expanding the amplitudes to the forth power of the momenta or pion masses 
(see [^ for an exception to this rule). Thus, the region of applicability of xPT is restricted to 
the very low-energy regime. In fact, xPT does not satisfy exactly the unitarity conditions but 
only in the perturbative sense, and this problem becomes very relevant at higher energies 0. 
To avoid this limitation of xPT several methods have been proposed like the introduction of 
new fields corresponding to resonances 0, the Pade approximants 0, the inverse amplitude 
method p, |^, etc. More recently it has also been considered the so called large N expansion of 
xPT (A^ being the number of Nambu-Goldstone bosons). This approximation is not restricted 
to the low energy domain and from the point of view of unitarity it has very good behavior. 
Until now this approach has been used for the computation of the pion scattering amplitudes 
[^ ^ in the chiral limit and also it was used in |n]| to study the 77 -^ tt'^tt" reaction. Some of 



the results in M were also reobtained in |11| but using different techniques. The work in |10 



has been, we think properly, criticised in |]T2[ , but work is in progress in order to make a correct 
computation of the 77 -^ tt'^tt^ reaction in the large A^ limit. 

The main goal of this work is to extend the work in p, |^ outside of the chiral limit i.e. we 
study the pion effective action and the pion scattering in the large A^ limit taken into account 



the effect of the pion mass. As it is well known this is not trivial at all since, because of the 
chiral symmetry, the introduction of the pion mass give rise to new interactions. 

The plain of the paper goes as follows: In section two we introduce the dynamics of the 
0{N + 1)/0{N) Non-Linear Sigma Model (NLSM) and we compute the effective action up to 
order 1/A^ including the leading corrections coming from having a pion mass different from zero. 
In section three we obtain the renormalized elastic pion scattering amplitude in terms of the 
renormalized coupling constants. In section four we study the running of those couplings and 
the renormalization of the pion mass and the pion decay constant. In section five we discuss the 
unitarity properties of the amplitudes. Section six is devoted to the particular case of the Linear 
Sigma Model (LSM) which is a very nice example to test the previously developped formalism. 
In section seven we show a fit of the experimental data for the I = J = elastic scattering 
amplitude and finally the conclussions are presented in section eight. 

2 The dynamics of the 0{N + 1)/0{N) = S^ non-linear 
sigma model 

In order to define the large N limit of xPT we start from the two-fiavor chiral symmetry group 
SU{2)fi X SU{2)l, and then the equivalence of the coset spaces SU(2)l x SU{2)r/SU{2)l+r = 
0(4)/0(3) = S^ is used to extend this symmetry pattern to 0{N + 1)/0{N) = S^ . Therefore 
A^ is the dimension of the coset space and the number of Nambu-Goldstone bosons (NGB), or 
in other words, the pions. The NGB fields can be chosen as arbitrary coordinates on S'^ . The 
most general 0{N + 1) invariant lagrangian can be written as a derivative expansion which is 
covariant with respect to both the space-time and the S^ coordinates. The lowest order is given 
by: 

^0 = ^^a6(7r)9^7r"9^7r^ (1) 



When using the standard coordinates the metrics is found to be: 

gab - dab + TTE^- ^^ l^J 

which has been obtained from the free lagrangian for the A^ + 1 fields vti, 7r2, ...ttat, a by forcing 
them to hve in the sphere vr^+cr^ = S^^VTaVTa+cr^ = NF^. As it is well known this lagrangian not 
only contains the kinetic term but also an infinite number of interacting terms with an arbitrary 
even number of pions. If we now want to introduce the pion mass we have to explicitly break 
the 0{N + 1) symmetry while keeping the 0{N) one. This can be achieved just considering in 
the lagrangian a properly normalized term proportional to the a component. Then the total 
lagrangian becomes: 



C = Igabd.n'^d'^n'' + NF'm'^l - ^ (3) 

where m plays the role of the pion mass as it can be seen by expanding the squared root. In 
addition we will also find another infinite set of terms which produce new pion interactions. 

Therefore the lagrangian in eq.3 will be our starting point to study the pion dynamics in 
the large A^ limit. As usual the action is given by: 

S[n]= f (TxC (tt) (4) 

However, to define the quantum theory we need to know not only the action but also the 
measure in the field space which defines the path integral representation of the Green functions. 
For the case of the NLSM, the appropriate measure includes a factor which is the squared root 



of the coset metrics determinant |T^. This factor contributes to the lagrangian with a term 



proportional to (5"(0). However, it is well known that using dimensional regularization |T5 



where n = 4 — e, it is valid the following rule, (5"(0) = or equivalently / (Pk = 0, and therefore 



the measure term can be ignored (for a recent discussion about regularization methods in xPT 
see 

Thus the generating functional for the regularized Green functions will be given by: 



d' 



^W= /[rf7r]e<'^W+/'^"^-^'"'^'") (5) 



with n = 4 — e. 

As it was discussed in the introduction we are interested in the computation of the effective 
action up to order 1/A^. Therefore we need some systematic method to obtain the diagrams 
that contribute in this approximation. This can be done using a generalization of the technique 
discussed in ||1^ for the case of the linear model, which introduces two auxiliary fields B^ and 
(f). The first one is connected with the derivative interactions whereas the second is related to 
those interactions proportional to m?. The lagrangian including these auxiliary fields is defined 
as: 



c = c-^-nfAb.. ^"'>'^' 



2 \ ^-FVi - im 



whose main properties are the following: first of all, since there is no kinetical term for the 
new auxiliary fields, they can be integrated out in an easy way, and then one finds: 



Second, the new self-interaction terms appearing in L have been chosen so that they cancel 
the old ones appearing in L and thus the new action can be written as follows: 



S[7i,B,<P] = f d^'x I NF^m^ --n'' (a + 771^)71"" 



^ AT 7-^5 / ri-J 2l2\ ^ ^^^^-^C:>o 7->// r I 



- -NF' [B' - m'(t)') - -NF'df.Bf'f 



7r2 



[nf^ 



V2NF^m^<Pg (-^] 



(8) 



where: 



fiv) = 2 l-Jl-77 



V 



9iv) = [^-i-^^-v) (9) 

This action describes the same dynamics for the pions than that in eq.3 and eq.4 but now, 
instead of self-interactions, we have B^ and mediated interactions. This fact makes it easier 
to trace out the A^ factors for a given diagram, and therefore to select those corresponding to 
the different terms in the Green functions 1/A^ expansion. 

From the action in eq.8 it is an straightforward matter to obtain the Feynman rules. To 
select the relevant diagrams for the leading order in the 1/A^ expansion the most important 
thing is to know the A^ power of the different propagators and vertex involved. It is immediate 
to check that the B and propagators are just constants proportional to 1/A^, the B2'k and 
(f)2iT vertices are of the order one, the BAn and the 047r vertices are of the order 1/A^, the BGir 
and the 067r vertices are of the order 1/A^^ and so on. 

The most general one-particle- irreducible {IP I) pion Green function can be obtained by 
attaching pions to the leg of the most general IP I B and Green function (see fig.l). Those 
IPI B and Green functions (represented by a black dot) can be expanded in terms of other 
reduced B and IPI functions (represented by a dashed dot) which are defined as the sum of 
all the diagrams containing only pion loops but not B or (p loops (see fig. 2 for an example). The 
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reduced B and Green functions are as much of order N (see fig. 3). In addition, the diagrams 
with B OT (f) loops contributing to the general B and (p Green functions are suppressed in the 
1/A^ expansion (because of the 1/N factors in the B and propagators) when compared with 
the first term. Thus we find the symbolic equation shown in fig.4. As usual, the connected 
Green functions can be expressed in terms of the IP/ ones by making trees, so that only the 
latter will be considered in detail. 

The B and IPI (dashed) Green functions can be obtained from the effective action: 



e^rr 



[^''^]= /[rfvr]e^/^"^^(-'^'^) (10) 



In fig. 5 we show the diagrams contributing at leading order in 1/A^ to those IPI Green functions 
with just one or two B and legs. It is not difficult to see that diagrams with three or more B 
or legs do not contribute at the level of precision considered in this paper, i.e. to the leading 
order in the 1/N expansion and to the order m'^/F'^. 

Now, using standard functional technics, it is possible to find: 

r[E,0] = frxlNF^m'--NF^[B^-m^(t)^)\ 

- V2NF^m^g(j^) Jrx(f){x) 

- 2Nm^g" (^) J dTxryc^ {x) K (x - y) <!> (y) 

- j/" (^) / rxd-yB, (x) G^-^ {x - y) B, (y) 

- tV2Nm'f' (^) g' (^) / rxd^yB, (x) G^ (x - y) <P (y) 



+ ... (11) 

where the points mean terms of 1/N^ order as well as terms with more than two B or (j) 
fields. Notice that the above effective action is not local being the integral kernels defined by: 



Gi,u {x - y) 
Gf^ (x - y) 

where for short we have defined 



/ ClrCC rvnrvi/1 I rZ 



dke-'^^'^-y^k^I (A;2 



(12) 



dk = fi^ 



d^~'k 



and the loop integrals /„ and /(A;^) are given by: 



(13) 



dq 



I (^e^ = I dq 



g2 _ ^ 



-m^ A 



{k + q) — rn^ {q^ — rn^) 
Using dimensional regularization these integrals are found to be: 



(14) 



and 



-m 



(47r) 



N, + I -log- 



in 



/i^ 



(15) 



lik 



2<|iV, + 2 + ^/l-^:;^log 



(47r 



A;2 



J- p J- , m 

log- 

+ 1 ^^' 



1 



Am? 
fc2 



(16) 



where, as usual: 



2 

iVe = - + log47r-7 (17) 

e 

Now we couple the effective action for the B and $ fields to the pions by defining the action: 



- V2NF'm'<j>g(-^]\+T[B,<j>] (18) 

As discussed before, the pion Green functions can be obtained by attaching pions to the B and 
$ Green functions. However it should be noted that it is also possible to have pion loops such 
as the one shown in fig.l, which are proportional to the integral /„ and therefore, to m^. Thus, 
from the formal point of view, the pion effective action can be written as: 

e'^^-^/zW = f ^dn'] [dB] [d0] e'-^'l"''^'^] (19) 

where the integral on the B and $ fields must be done at tree level. Moreover, since we want 
to expand the pion effective action only up to order m'^/F'^ terms, we only have to perform the 
integral on the pion fields to the one loop level. Therefore, to the considered degree of accuracy, 
we have: 



5[7r',5,0] = S 



7l',B,i 



+ ^/'""'"^ ^.-(xwm '"-"^"^"^"^^^^^^---- ^'°^ 



[x) Stt' (y) 
where B, and vr are defined so that: 



ss 

SB^ {x) 

SS 

6(f) (x) 

6S 

Sn^ (x) 









(21) 



for B = B, 
follows: 



and 7t' = 71. These equations can be used to write B and in terms of tt as 



and 



BA^) 



(fydke 



-ik{x—y) 



ik,^ 



i + ^p^ 



n 



X 



(y) 



I (k) rn^ 



2NF^ 



1 + 



k^I{k) 
2F2 



F2 



l + A 



F2 



[X 



dTydke 



-ik{x—y) 



1 



iVF2 ^ 4 V NF^ 



'^HvY 



k^I{k) 
2F2 



1 



k^ijky 

2F2 



+ 



(22) 



(23) 



If we now integrate the B and $ fields to the tree level and the pion field to the one-loop 
level we can write: 



S, 



eff 



TT 



s 



vr,5,0j+-Trlog^-;,^-^^-^ 



+ 



(24) 



where: 



10 



,-,1^ . .ia,^-r^S.,i.-,n20^Br{l^y-^^ 



+ d,B^f'[^^j5at6ix-y) + ... 
= Dab{x,y) + Aabix,y) 



(25) 



with 



Dab{x,y) = -(n + m^)6abS{x-y) (26) 

and Aab{x,y) can be obtained from the above two equations. Now it is possible to write: 



Tr log (D + A) = Tr\ogD[l + D-'A 



°° f-l 
Tr log JJ + y^ 



,«+! 



n=l 



n 



D^^A 



(27) 



In the following we will ignore the Tr log D term since it does not depend on the pion fields 
and therefore its contribution to the pion effective action is irrelevant. On the other hand it is 
not difficult to see that only the term n = 1 contributes to the effective action up to the level of 
the approximations considered here. Then, using the above equations it is possible to arrive to: 



^eff [vr] 



+ 



X 






SNF'^ 
1 



ct'xct^y 



(27r) 
m2(l + A|^) 



^'^("^^^2 {x) TT^ {y) 



1 + 



+ O 






O 



1 + 



m 



2F2 



(28) 
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which is our final result for the regularized pion effective action to leading order in the 1/A^ 
expansion and to order m? / F"^. 

3 The Scattering Amplitude 

immediate to find the scattering amplitude for the process liaT^b -^ T^c^d from the pion effective 
action obtained in the preceding section. This amplitude can be written as: 

Tahcd = Sab5cdA (S) + SaAdA (t) + 6adSbcA (u) (29) 

where: 



A(s) 



NF"^ 



w? ( ^ . s 



1 + ^ fl + £^fV F2 



1 + A™ (30) 



2F2 

It is not difficult to see that the above amplitude corresponds to the diagrams appearing 
in fig. 6, which are those relevant to the leading order in the 1/A^ expansion and to order 
m? / F"^. However, for practical applications, this regularized amplitude should be renormalized. 
Moreover, the renormalization procedure is not trivial since the NLSM is not renormalizable in 
the standard sense. As it was discussed in the renormalization of the NLSM in the large 
N limit requires the introduction of an infinite number of counter terms in the renormalized 
lagrangian. In the case considered here things are even harder since, in addition to the derivative 
interactions that appear in the chiral limit, we also have new interactions proportional to iv?. 

Therefore, we will introduce the renormalized lagrangian: 

Cr = C + c.t (31) 
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where c.t. represents the infinite set of counterterms needed to absorb all the divergences 
appearing in the the effective action in eq.28. In fact, at the level considered here, we do 
not need to know the precise form of these counterterms. From fig. 6 we know that the only 
vertices appearing in the diagrams contributing to the regularized pion scattering amplitude 
are the derivative four pion vertex, the four pion vertex proportional to m^ and the derivative 
six pion vertex. The effect of the counterterms needed for the renormalization of the pion 
scattering amplitude on these vertices can be easily taken into account just making the following 
replacements: 



q^ q^ 



on the derivative four pion vertex 



G (32) 



2 2 

m m 



NF"^ NF"^ 

on the four pion vertex proportional to m^ and 



H (33) 



(34) 



{NF^Y {NF^ 

on the derivative six pion vertex (in the above prescriptions q represents the total momentum 
of two of the pions involved in the vertex although for the sake of simplicity we do not give 
the details of the appropriate combinatory ). The factors G, H and E are arbitrary analytical 
functions on s so that they can be written as: 



s f s \'^ °° /' s ^'^ 

G{s) = ^+9i-E^ + 92(-e;^) +... = Y.9k 



k=l 
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His) = i+h4,+hJ4,)\... = f:hJ4,Y 



E{s) = l + ei^ + e,(^)V... = |:e.(^)' (35) 

The constants gk, h^ and Ck are the unrenormahzed couphng constants and must be renor- 
mahzed to absorb the divergencies in the pion scattering amphtude. Here it is important to 
notice that, as it can be seen in eq.28, the regularized action has not any new contribution to 
the kinetic term, and therefore no renormahzation of the pion wave function is needed. Now 
we can use the above replacements to compute again the diagrams in fig. 6. The new amplitude 
has the same form that the previous one with 



1 


sG{s) 


2 

m 


H (s) + AE (s) -^' 


NF^ 


1 1 sG(s)I{s) 




r. 1 sG{s)Iis)] 


2 



A{s) = 

= Aois) + A,{s) (36) 

where Ai{s) is the term proportional to m^ and Ao{s) = A{s) — Ai{s). As there is no 
renormahzation of the pion wave function, the above amplitude should be made finite by an 
appropriate definition of the renormalized couplings in terms of the unrenormahzed ones. This 
can be done as follows; first we will consider the renormahzation of Aq and Ai separately. For 
the renormahzation of the Aq amplitude we follow similar steps to those followed in p[ for the 
renormahzation of the pion scattering amplitude in the chiral limit. Thus we write: 

NF^ N N 

where from eq.l6 and eq.l7 T{s; fi) is given by 



4777.2 , V 1 T 1 . m'^ 



r(,;^)^,/l-— logJL=^_-,og- (38) 
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Now we define G^ as: 



(iV, + 2) 



GR G 2(47r)'F2 
The G^ function can be expanded in terms of the renormahzed coupling constants gj}: 



(39) 



G« = l+gf(^)+ft«(^)%... 



s 
F2 



(40) 



so that, expanding eq.39 in powers of s/ F"^ we can find an infinite number of equations 
defining the renormahzed couphngs g^ in terms of the unrenormahzed ones g^- 
With the above definition for G-^, we can write: 

1 sG^ 



An(s] 



NF^l 'GR T(s-u) 



(41) 



2(47r)^F2' 

In this form Aq is a manifestly finite and well defined function on s once G^ is given. 

The next step is the renormalization of the second piece of the amplitude Ai. First we note 
that by comparison of eq.36 and eq.41 we can write 



so that Ai is given by: 



G' 



Alls) 



G 






-m'(^) (h + AE^ 



(42) 



2F2 



T{s;fi) 



(43) 



Taking into account that the Aq in eq.41 is finite we also have to require: 



1 
G2 



H'- 



Es m 

/i2 



{A-kY F2 



(44) 
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to be finite, where: 



Es 



H' = H+j—^iN, + l) (45) 



In other words, the quotients H'/G^ and E/G'^ must be finite. Therefore we introduce the 
two new generating functions A^ and B^ follows: 

H' = A^'G^ 
E = B^G'^ (46) 

which generate two infinite set of renormalized constants: 

A« = l + af-^ + ... 

B^ = l + bf^^+... (47) 

Thus, the A{s) function appearing in the elastic scattering amplitude can be written in terms 
of functions which are finite in the e ^ limit. 



Ais) 



1 




sG^ 


m^G''^ 


[A"" - B^j^hgf] 


NF^ 


1- 


-I0k^ns;f^) 




1 


-I0k^ns;f^)] 


2 



(48) 



4 Renormalization Group Equations 

In the last section we have been able to write the scattering amplitude in a manifest finite form 
in terms of the G^, A^ and B^ functions. These functions generate three countable infinite sets 
of renormalized couplings which implicitly depend on the renormalization scale fi. Therefore, 
G^, A^ and B^ also depend on /i. As we have already noticed, there is no wave function 
renormalization in the leading order of the 1/N expansion, and thus we conclude that the A{s) 
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amplitude is a physical observable in the sense of the renormalization group evolution. In other 
words, the explicit dependence of A{s) on the renormalization scale n must be exactly canceled 
by the implicit dependence through the generating functions. This can be precisely stated with 
the evolution equation: 



Q I 3 3 3 



A{s) = (49) 



where, as usual, the beta functions are the derivatives of the corresponding couplings with 
respect to log/x. The above equation is extremely useful since it can be used to determine the 
dependence of the generating functions on /i, and finally, the dependence of all the couplings 
they generate on this scale. In fact, this can be done in many different ways but, in order to 
make the chiral limit more transparent, we will require the above equation to apply separately 
to Aq and Ai. Thus, from eq.41 and eq.49 we can write 



dG"" (s; fi) (l - £^i!i£^T (s; fi)] + G^ (s; fi)d(l- f^^^ifi^T (s; fi)] = (50) 



2(47r)"F2 ' '"7 ' '"' V 2(47r)'F2 

and integrating this equation we find: 



2{47r)^F2 ^& MO 

This equation describes the dependence on the G^{s; fi) generating function on the renor- 
malization scale /i determined by the renormalization group equations. 
Now, applying eq.49 to Ai we find that the piece 



2 

A^{s- /i) - B%s- /i) . ^ log — (52) 
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must be n independent. For the sake of simplicity we introduce another generating function 
J-^(s;/i) defined as: 



2 2 

J^(s; a) = A^(s; a) - B^(s; a) % log ^ + % log ^ (53) 

note that in general the couplings jj} generated by J^ can also depend on logm^ but this is 
not the case of jf . In terms of J^ eq.52 reads: 



dJ^{s;i2) _ -2s 
d{logfi) ~ (47r)^F2 



(54) 



This equation can be integrated to give: 



J^ (s; /i) = J« (s; /io) - 77^^ log - (55) 

[4:71) F^ /io 

Now it is possible to expand both sides of the evolution equations eq.51 and eq.55 for the 

G^ and the J^ generating functions in powers of s/F"^: 

G"" {s; f.) = Y. 9^ if^) (i^) (56) 

fc=o ^-^ ^ 

fc=0 ^-^ ^ 

to find two infinite sets of evolution equations which explicitly give the dependence of the 

renormalized couplings generated by G^ and J^ on the scale /i. For example, for the lowest 

couplings we find: 



(7f(/i) = (?f(/io)--— logi^ 



jf (-") = Jf (/io) - A log ^ (57) 



w" /io 
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In terms of the G^ and J^ generating functions the A{s) amphtude can be written as: 



A(s) 






2{47r)^F2 - 



J\s-fi)- 



o log- 

{Anf F^ fi 



(5^ 



From this result it is not difficult to obtain the the amphtude low energy behavior which is 
given by: 

i2 



A is) 



NF^ 
O 



1 



F2 J 






■m 

iVF2 



(59) 



This is a very useful and interesting result since we know that at low energies A{s) should go as 
(s — m\)/ f"^ where 771^^ is the physical pion mass and /^ is the pion decay constant. Therefore 
the amphtude in eq.58 has the right low energy behavior provided we define: 









NF^ I 1 - 



m 



,R 



1 TTl 



,2 



R 



o 



m 
F2 



(60) 



this equation also provides the right non trivial dependence of f^^ on the logarithm of the pion 
mass which is well known from chiral perturbation theory and in it F'^pf is not depending on 
the scale /i as can be shown by using eq.57. These facts are therefore a good consistency check 
of our results. Moreover, from eq.28 we see that the pion mass should also be renormalized. 
This can be done, for example, defining the renormalized mass as: 

,2 "1 



m 



R 



m^U 



m 



{N. + 1) 



2(47r)"F2 
so that the effective mass appearing in the effective action can be written as: 



(61) 



m^ 



m\ + 



m 



, m 
log^ 



(62) 



'f^ "^^ ' 2(47r)^F2"">2 
in good agreement with well known results from chiral perturbation theory ||^. However it is 
important to notice that our computation of the effective action and the scattering amplitude 
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for pions is only performed up to the lowest order in m? /F"^. Thus, as the renormalization 
of the mass introduces an extra factor m? / F'^ i^eff — ^ + O (^j), we do not have to 
renormalize that m? appearing in the amplitude in eq.30 since it would have produced extra 
terms proportional to {m? / F'^Y, and the same applies to eq.58. Therefore we can consider the 
rn? appearing in eq.30 and eq.58 as the physical pion mass squared m^. 

5 Partial Waves and Unitarity 

In order to study the unitarity properties of the amplitudes obtained in the preceeding section 
we will perform the standard partial wave expansion. First we project on the isospin channels 



which for the model here considered are be defined as [|TT|: 



To = NA{s) + A{t) + A{u) 

Ti = A{t)-A{u) 

T2 = Ait) + A{u) (63) 

The partial waves are then defined as usual: 

ajj(s) = [ d (cos 9) Ti is, cos 6) Pj (cos 6) (64) 

647r J-i 

As it is well known, the requirement of unitarity constraints the possible behaviour of these 
partial waves. In particular, they should have a cut along the positive real axis from the 
threshold to infinity. The physical amplitudes are obtained from the values right on the cut, 
and, in addition, the condition of elastic unitarity: 

Imajj = cr|a/j|^ (65) 
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(where a = J\ — Am? / E'^) must be fulfilled in the physical region s = E'^ + ie and E'^ > Aw? 
where E is the center of mass energy. This equation is exact for energies below the next four 
pion threshold but even beyond that point it is approximately valid. 

In the following we will study how well this relation is satisfied in the leading large N pion 
scattering amplitudes. From eq.63 it is obvious that the most important isospin channel at the 
leading order in the large N expansion is / = which is of order 1 in the 1/A^ expansion. From 
eq.58, eq.63 and eq.64 and for the aoo channel it is immediate to find: 

1 I sGr{s-ii) ^'('/i?(s;/i)-(i^iog^] 



aoo {s} 



327rF2 1 _ sGnjsp,) (f ■ 



fiS^r^+-- 



#(t + 



(66) 



where we have defined T as 



T [E +ie;fj, ] = fflog -\+i7ia -log^ 

= f^mo (67) 

and thus we are only considering explicitly the physical values of s. However, it is easy to see 
that the aoo(s) partial wave amplitude, when s is considered a complex variable, has the proper 
right unitary cut mentioned above. The partial wave above can be now naturally decomposed 
as: 

^00 = ctoo + '^00 (68) 

where Cqo corresponds to the part of the partial wave that is proportional to vr? and agQ is the 
part which is not. In terms of ago and ago the unitarity condition in eq.65 reads: 

Jmago = o"|aooX 

ImaQQ = 2(Ti?eaoo Ooo (69) 

provided one neglects terms of order {ir? / F'^Y. After a bit of algebra it is not difficult to show 
that these equations are indeed satisfied by ago and a^Q as obtained from eq.65 and eq.68. 
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Therefore, we finally can write for the partial wave aoo from the amplitude in eq.66 : 



ImaoQ = o-|aooP + 0{ — ]+0 



m 
F2 



(70) 



As usual the right cut produces two sheets for the aoo function. Depending on the actual form 
of the generating functions, eventually some poles could appear in different places. If the poles 
appear in the first sheet (the physical sheet) they are not acceptable and must be understood 
as artifacts (ghosts) of the approximation. However, when the poles appear on the second sheet 
(the unphysical sheet) they are welcome and have a natural interpretation as resonances. This 
is for example the case of the LSM (Linear Sigma Model) that will be studied in detail in the 
next section. 

6 The Linear cr-model 

In this section we consider the case of the LSM which is defined by the lagrangian 

C = ^d^(l)''d^(l)-V{(l)) + CsB (71) 

where is an A^ + 1 vector with components 7ri,7r2...,7r7v,7rAr+i and the potential is given by: 

\/(0) = VI0P + AI0T (72) 

The lagrangian would be 0{N + 1) invariant were it not for the last symmetry breaking piece 
that it is only 0{N) invariant. This piece is given by: 

Cgs = ^/NFm^a (73) 

where 

a = TT^v+i = ViVF2 + H (74) 
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In terms of the field H the lagrangian can be written as: 

2 M 2 ^ 

- AXNF^H^ + HVNF^m^ -X(n^ + H^^ ^ 



- AXVNF^Hti'^ -AXVNF^H^ (75) 

Notice that there is a hnear term in the H field. In order to eliminate it, we introduce a shifted 
h field: 

H = Ho + h (76) 

with constant Hq. From the lagrangian written using this new field h, and in particular, from 
the term proportional to tt^ we find the pion mass to be: 



m 



2 = 4XH^ + SXHoVNF^ (77) 



The piece of the lagrangian proportional to h can be written in terms of this pion mass. Then, 
the condition for this piece to vanish become: 



m'^ 



Ho = -VNF^ + \NF^ + — - > (78) 

V 2,X 



So that the lagrangian has the final form: 



2 M 2 2 ^ 2 



,2 . , , / o , o\ / ,^„„ Vfl^ 



X (tt^ + /^y - 4A/1 (vr^ + h") JnF^ + ^ (79) 



where 



ml = M^ + 3m^ 

M^ = 8XNF^ (80) 
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Now we are able to obtain the corresponding Feynman rules. The elastic scattering amplitude 
can be computed in many different ways in the large N limit. Probably the easier one is the 
following: First one compute the diagrams in fig. 7a, then one iterates the result as shown in 
fig. 7b and, finally, one expands in powers of m^ / F'^ and retains the zero and the first order. 
The final result is given by: 

sGm{s) m2Gi,(s)(l + 2-^ + A^^ 



j\_(g\ 



where: 



Gm (s) = -— ^ (82) 

Notice that we can consider this result is a particular case of the general one given in eq.36 with 
the elections: 

G{s) = GMis) 

His) = Gi,(s)(l + 2^) 

E{s) = Glis) (83) 

The possibility of writing the LSM amplitude as a particular case of the NLSM (in the large A^ 
limit in both cases) was expected form our previous discussion but, as far as both computations 
were done in a completely different and independent way, it turns to be a highly non trivial 
check of our results. Now we can proceed with the renormalization procedure as we did in sec. 3 
and sec. 4. The relevant issue here is that the LSM is renormalizable in the standard way, i.e. 
only a finite number parameters have to be renormalized. This means in the large A^ limit here 
considered that only the constant M will be renormalized. In particular, we find: 

G^ {s; li) = \ (84) 

M^(s;m) 
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where 

1 1 



ms-m= + J(3;jV^W + 2' («=) 



so that the dependence of the M^ with the scale is given by 



The other needed renormahzed generating function is: 



MUi^) = ^^ (86) 

■■- 9.(dZ\2p2 log -3 



^''(^^''> = ^ + 2m|W '''' 



Now, using eq.86 we find: 



dJji (s; /i) —2s 



t/(log/i) (47r)V2 ^ ^ 

which is consistent with the general /i dependence for the generating functions found in eq.51 
and eq.55. Finally, the renormahzed amplitude for the LSM is found to be: 



2(47r)^F2 V '/-; (^1 2(47r)^F2 



In the chiral limit, this amplitude reduces to that found in |]T8[. The I = J = partial 
wave amplitude has a pole in the second Riemann sheet that corresponds to a physical scalar 
resonance. This resonance has the interesting saturation property, which means that its mass 



remains bounded when Mji goes to infinity (see also |T9[ for a general discussion). This fact is 



probably related with the expected triviality of the LSM. 

In any case, the LSM provides a very nice example of how the general formalism developed 
for the computation of the large N limit of xPT works outside of the chiral limit. 
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7 Fitting the pion scattering 

Apart from the LSM, the most interesting case to apply our results is the pion scattering. 
However, in order to use eq.58 to fit the experimental pion scattering data, we have to face 
the problem of dealing with an infinite number of parameters, i.e. the scale fi and the values 
of the renormalized coupling constants g^{fi) and Jk{fi) at that scale. However one can easily 
solve this problem just considering only particular cases where all the coupling constants but a 
finite set g^, g2 ■ ■ ■ g^ and jf , J2 ■ ■ ■ j^ vanish at some scale fi. These models are just defined 
by a finite number of parameters (/i and the r + s coupling constants renormalized at this 
scale) and therefore can be used to fit the experimental data. In particular, one can consider 
the extreme case of having all the renormalized couplings equal to zero at some scale /x, i.e. 
9k if^) — Jkif^) — f*^^ all fc > (the minimal model). The model thus obtained only has one 
parameter (the /i scale) and it can be considered in some sense as the non-linear sigma model 
renormalized at the leading order of the 1/A^ expansion out of the chiral limit. 

In fig. 8 we show the result of our fit of the I = J = phase shift for elastic pion scattering 
done with the above defined minimal model. The only parameter in this specially simple case 
is the scale /i. The fit in fig. 8 corresponds to a value of /i ^ 775MeV and as it can be seen it 
describe perfectly well the data from the threshold to 700Mev. It is also interesting to remark 
that the fitted value for fi is the same that was found in |^ . However, as we are including here 
the effect of the pion mass, the fit is much better in the threshold region. 

For other channels like / = J = lor/ = 2, J = the prediction of the leading order of the 
large A^ expansion is that they are supressed. This could seem quite a poor prediction but it 
is not. Looking at the experimental data we can realize that the phase shifts in these channels 
grow very slowly with energy, for instance, at 500MeV of center of mass energy we have ^n ~ 5 
degrees and ^20 — —7 degrees whereas 5oo — 38 degrees at the same energy. In any case one can 
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use the results obtained in the previous section considering eq.58, eq.63 and eq.64 with the /x 
value fitted for the aoo channels. The results are shown in fig.8 and fig. 9. For the / = 2, J = 
channel we see that our fit (wich is in fact a prediction since we only have one parameter that 
had already been fixed in the I = J = channel) is quite good. However this is not the case 
when I = J = 1. Basically this channel can be described as a non interacting channel at low 
energies but a strong interacting channel at higher energies due to the appearence of the p 
resonence. The approach considered here seems to describe well the non-interacting low-energy 
region but it fails in the resonant region. One is then tempted to make an interpretation of 
the scale fi as some kind of cutoff signaling the range of applicability of our approach and, in 
fact, this interpretation had already been discussed in p. From this point of view, it is quite 
interesting to realize how close the /i fitted value (/i ~ 775MeV) is to the p mass. In some 
way, by fitting the I = J = channel, the model is telling us where new physics can appear in 
other channels like the I = J = 1 (the p resonance) and thus setting the limits of applicability 
of the model. However, we would like to stress that our renormalization method is completely 
consistent and our results are formally valid at any energy independently of the goodness of the 
fits we can obtain with them. 

Finally, for comparison we have also included in fig.8 and fig. 9 the results obtained with 
standard xPT to one loop with the parameters proposed in [^. It is important to note that, 
in spite of using many more couplings, it does not a better job than the approach considered 
here. The description of the I = J = channel is good in both cases, the J = 2, J = is much 
better in the large A^ limit and the J = J = 1 is good at low energies (where the channel is not 
interacting) but it fails to describe the resonance region too. 
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8 Conclusions 

In this work we have study how the large N approximation (A^ being the number of Goldstone 
bosons) can be apphed to xPT in order to improve its unitarity behavior and enlarge its energy 
applicability region. In particular we have computed the renormalized effective action at the 
leading order including the effects due to the pion mass up to order in\/ f"^. The amplitude 
obtained from this action has the appropriate low energy limit and provides the right dependence 
of the the pion mass and the pion decay constant on the chiral logs. We use this action to 
reproduce the LSM and to fit the vr — tt scattering with just one parameter. Within this 
approach the fit is competitive with the standard one-loop fit both failing in the description of 
the p resonance. However, it can improve when next to leading (1/A^^) corrections are included. 
Note that due to the structure of the large N amplitudes it is possible to find poles in the second 
Riemann sheet that could reproduce resonances as it happens in the large N description of the 
LSM. At present, work is also in progress in that direction ||13|| . 

In conclusion, we consider that the great sucsess of xPT can be enlarged with the use of 
well defined non-perturbative techniques as the large A^ expansion. The improvement can be 
obtained in several directions like unitarity, energy range of application and finally, in order to 
make predictions, avoiding the inflation in the number of free parameters. 
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Figure Captions 

Figure 1: One typical Green function obtained by attaching pions (continous line) 
to one B (wavyline) and $ (dashed line) Green function. Note the possibility of 
having pion loops in the attaching point. 

Figure 2: Example of $ Green (black dot) in terms of the reduced (i.e. including 
only pion loops) $ and B Green functions (dashed dot). 

Figure 3: Example of leading diagramns contributting to the reduced B Green 
functions in the Large N limit. 

Figure 4: Relation between the B Green functions and the reduced B Green 
functions at leading order in the large A^ approximation. 

Figure 5: a) Diagrams contributing to the reduced one $ Green function b) Di- 
agrams contributing to the reduced two B Green function. Similar diagrams con- 
tribute to the two $ and the 5$ reduced Green functions. 

Figure 6: Diagrams contributing to the the pion scattering to leading order in the 
large A^ approximation up to order im? / F'^. The large black dot represents the pion 
coupling proportional to m? 

Figure 7: Diagrams appearing in the computation of the pion scattering amplitude 
in the LSM described in the text. The continuous lines represent pions and the 
dashed ones the h field. 

Figure 8: Phase shift for vrvr scattering. The results coming from the large A^ limit 
out of the chiral limit are drawn with a dashed line both for the (0, 0) and the (2, 0) 
channels. The continuous lines were taken from [^ and they are the (0, 0) and (2, 0) 
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phase shifts as obtained from the standard one-loop xPT with the parameters given 
in 0. The experimental data corresponds to: [^ A, [g2[ Q, [^] D, |2g ^, |2§ 
X, and ||28ll •• 



V, HI*, 

Figure 9: (1, 1) Phase shift for tttt scattering. The dashed line corresponds to our 
fit using the large A^ limit. The continuous line |^ is the result coming from one-loop 
XPT with the parameters proposed in [^ . The experimental data comes from: 

o, m A. 
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